Magnetic-field induced spin-Peierls instability in strongly frustrated quantum spin 

lattices 
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For a class of frustrated antiferromagnetic spin lattices (in particular, the square-kagome and 
kagome lattices) we discuss the impact of recently discovered exact eigenstates on the stability of 
the lattice against distortions. These eigenstates consist of independent localized magnons em- 
bedded in a ferromagnetic environment and become ground states in high magnetic fields. For 
appropriate lattice distortions fitting to the structure of the localized magnons the lowering of mag- 
netic energy can be calculated exactly and is proportional to the displacement of atoms leading to a 
spin-Peierls lattice instability. Since these localized states are present only for high magnetic fields, 
this instability might be driven by magnetic field. The hysteresis of the spin-Peierls transition is 
also discussed. 

PACS numbers: 75.10.Jm, 75.45.+j 
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Antiferromagnetically interacting spin-^ systems on 
geometrically frustrated lattices have attracted much 
attention during last years. Such systems have rich 
phase diagrams exhibiting a number of unusual quan- 
tum phasesi*^. A striking example is the kagome lattice 
antiferromagnet having a liquid-like ground state with a 
gap for magnetic excitations and a huge number of singlet 
states below the first triplet state (see e.g. Ref. 13 and ref- 
erences therein). Another intriguing example is the struc- 
tural phase transition in spin systems driven by magne- 
toelastic coupling (spin-Peierls instability) observed, e.g., 
in CuGeO^. Frustrating interactions may also provide 
a route to generating fractional phases in two dimensions 
which manifest itself in the dynamic correlations probed 
by inelastic neutron scattering experimental. 

In the presence of an external magnetic field frustrated 
quantum spin systems exhibit a number of unusual prop- 
erties. In particular, plateaus and jumps can be observed 
in the zero-temperature magnetization curve for such 
modelsiSi. The theoretical investigation of exotic magne- 
tization curves has been additionally stimulated by the 
experimental observation of plateaus e.g. in CsCuCla^ or 
SrCu2(B03)2^. 

Though the treatment of quantum spin systems often 
becomes more complicated if frustration is present, in 
some exceptional cases frustration is crucial to find sim- 
ple ground states of product form^-^*'. Recently, for a 
wide class of frustrated spin lattices exact eigenstates 
consisting of independent localized magnons in a ferro- 
magnetic environment have been foundii. They may be- 
come ground states if a strong magnetic field is applied 
and lead to a macroscopic jump in the zero-temperature 
magnetization curve just below saturation. 
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In the present paper we examine the stability of an- 
tiferromagnetic spin lattices hosting independent local- 
ized magnons with respect to lattice distortions through 
a magnetoelastic mechanism. We are able to present rig- 
orous analytical results completed by large-scale exact 
diagonalization data for lattices up to = 54 sites. We 
discuss the field-tuned changes of the ground-state prop- 
erties owing to a coupling between spin and lattice de- 
grees of freedom. 

The efi^ect of a magnetoelastic coupling in frustrated 
antiferromagnets is currently widely discussed. First of 
all the investigation of the Peierls phenomena in frus- 
trated ID spin systems is in the focus (see, for instance, 
Ref. and references therein) . But also in 2D and 3D 
quantum spin systems lattice instabilities breaking the 
translational symmetry are reported. The frustration- 
driven structural distortions in the spin-| square-lattice 
J1-J2 Heisenberg antiferromagnet studied in Ref. ITsI 
might be relevant for Li2VOSi04 and VOMoO^iiii. 
For 3D frustrated antiferromagnets containing corner- 
sharing tetrahedra one discusses several examples for lat- 
tice distortions. Inelastic magnetic neutron scattering on 
ZnCr204 revealed that a lattice distortion can lower the 
energy driving the spin system into an ordered phasei^. 
NMR investigation of the three-dimensional pyrochlore 
antiferromagnet Y2M02O7 gives evidence for discrete lat- 
tice distortions which reduce the energjii^. A lifting of a 
macroscopic ground-state degeneracy of frustrated mag- 
nets through a coupling between spin and lattice degrees 
of freedom in pyrochlore antiferromagnets was studied in 
Refs. mil which may have relevance to some antifer- 
romagnetic compounds with pyrochlore structurei^. 

In all those studies the lattice instability is discussed at 
zero field. As pointed out already in the late seventies^S a 
magnetic field may act against the spin-Peierls transition 
and might favor a uniform or incommensurate phase. In 
contrast to those findings, in the present paper we discuss 
magnetic systems for which the magnetic field is essential 
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FIG. 1: Square- kagome lattice with one distorted square (left) 
and kagome lattice with one distorted hexagon (right) which 
can host localized magnons. The parts of the lattices before 
distortions are shown by dashed lines. All bonds in the lattices 
before distortions have the same length. 



for the occurrence of the lattice instability. 

To be specific, we consider two geometrically frustrated 
lattices, namely, the square-kagome lattice (Fig. ^ left) 
and the kagome lattice (Fig. right) . The ground state 
and low-temperature thermodynamics for the Heisenberg 
antiferromagnet on both lattices are subjects of intensive 
discussions (see e.g. Refs. fillll2Ci2l() . The Hamiltonian 
of N quantum {s ~ i) spins reads 

i^ = ^ Jnm ( ^ + S-S+) + As^^S^j ~ hS\ (1) 

(nm) ^ 

Here the sum runs over the bonds (edges) which con- 
nect the sites (vertices) occupied by spins for the lattice 
under consideration, Jnm > are the antiferromagnetic 
exchange constants between the sites n and m, A > 
is the anisotropy parameter (in most cases throughout 
this paper A = 1), /i is the external magnetic field and 
= J2n is z-component of the total spin. We 
assume that all bonds in the lattice without distortion 
have the same length and hence all exchange constants 
have the same value J. 

From Ref. ^3 we know that independent localized one- 
magnon states embedded in a ferromagnetic background 
are exact eigenstates of the Hamiltonian ^ for the con- 
sidered models. More specifically, by direct computation 
one can checkii that 

|1) = ^ E (--^Y^ilFM): (2) 

i^square 

|i) = ^ E (-i)x-i^A^) (3) 

i^hexagon 

are the one-magnon eigenstates of the Hamiltonian 
on the square-kagome and kagome lattices, where \FM) 
stands for the embedding fully polarized ferromagnetic 
environment. The corresponding energies {h = 0) of the 



one-magnon states (0) and ^ are 

- J+ J+(27V-12)ij, (4) 

-ij + 2J+(27V-18)ij. (5) 

The magnons ((2Jl or ||2J) are trapped (localized) on a 
square or on a hexagon, respectively. We separate explic- 
itly in Q , the contributions to the energy from those 
bonds which form a magnon trapping cell (first terms), 
from the bonds connecting this cell with the environment 
(second terms) and from the ferromagnetic environment 
(third terms). One can proceed to fill the lattices with 
n > 1 localized magnons; the state with maximum filling 
("magnon crystal") has n = n,„ax independent localized 
magnons with rimax — ^'Hd for the square-kagome 
and kagome lattices, respectiveljiii. Since commutes 
with the Hamiltonian (Q, the energy in the presence of 
an external field h ^ 0, E{S^,h), can be obtained from 
the energy without field, E{S^), through the relation 
E{S^, h) = E{S^) — hS^. Since each magnon carries one 
down spin, a localized magnon state with n independent 
magnons has S"^ = — n and consequently Eqs. I@J, 
© give E{S^ — — 1) for the corresponding systems. 

Under quite general assumptions it was proved^^ that 
these localized magnon states have lowest energies in 
the corresponding sector of total S^. As a result, these 
states become ground states in a strong magnetic field. 
More specifically, the ground-state energy in the pres- 
ence of a field is given by Eq^H) — Ej^miS^) — hS^ and 
the ground-state magnetization is determined from 
the equation h = E^in{S^) — E^in{S^ — 1). Since for 
S"^ = iiV, . . . , iiV - ninax (i.e., < n < Umax) the local- 
ized magnon states are the lowest states, one has 

E,^in{S^)^^NJ-3nJ = -NJ + 3JS' (6) 

for both models. Due to the linear relation between 
Sniin and one has a complete degeneracy of all lo- 
calized magnon states at h = hi, i.e. the energy is 
-NJ at h = hi for all ^A^ - n^ax < S"' < ^N, where 
hi = 3J is the saturation field (identical for both mod- 
els). Consequently, the zero-temperature magnetization 
jumps between the saturation value and the value 
(jgN) for the square-kagome (kagome) lattice. The 
effects of the localized magnon states become irrelevant 
if the spins become classical (s — > oo). 

We want to check the lattice stability of the consid- 
ered systems with respect to a spin-Peierls mechanism. 
For this purpose we assume a small lattice deformation 
which preserves the symmetry of the cell which hosts the 
localized magnon (in this case the independent localized 
magnon states remain the exact eigenstates) and analyze 
the change in the total energy (which consists of the mag- 
netic and elastic parts) to reveal whether the deforma- 
tion is favorable or not. To find a favorable deformation 
one needs optimal gain in magnetic energy. For that we 
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use the circumstance that due to the locahzed nature of 
the magnons we have an inhomogeneous distribution of 
nearest- neighbor (NN) spin-spin correlations (siSj)Ai. In 
case that one magnon is distributed uniformly over the 
lattice the deviation of the NN correlation from the ferro- 
magnetic value, (siSj) — i, is of the order -i.. On the other 
hand for a localized magnon |j2Jl [0] we have along the 
hexagon [square] hosting the localized magnon actually a 
negative NN correlation (siSj) — ^^nd all other 

correlations are positive. Hence a deformation with opti- 
mal gain in magnetic energy shall lead to an increase of 
antiferromagnetic bonds on the square (hexagon) and to 
a decrease of the bonds on the attaching triangles. The 
corresponding deformations are shown in Fig. ^ For 
the square-kagome lattice (Fig. ^ left) the deformations 
lead to the following changes in the exchange interac- 
tions: J (l + \/2S) J (along the edges of the square) 
and J ^ (l — (-\/3 — l) S) J (along the two edges of 
the triangles attached to the square) , where the quantity 
d is proportional to the displacement of the atoms and the 
change in the exchange integrals due lattice distortions is 
taken into account in first order in S. For the kagome lat- 
tice (Fig. m right) one has J ^ {I + S) J (along the edges 
of the hexagon) and J ^ (^1 — ^S) J (along the two edges 
of the triangles attached to the hexagon). The magnetic 
energies (0J and ((SJl are lowered by distortions and be- 
come ^NJ-3J- \V2 {3 + VS) SJ a.ud ^NJ-SJ-^SJ, 
respectively. This is in competition with the increase of 
the elastic energy, which is given in harmonic approxi- 
mation by 2 (6 — -s/S) ad^ (square-kagome) and by 9aS^ 
(kagome). The parameter a is proportional to the elastic 
constant of the lattice. The change of total energies due 
to distortions read 

-^V2(3 + VtjSJ + 2(6-VtjaS^, (7) 

^^SJ + 9aS^, (8) 

where for n independent localized magnons with n re- 
lated distortions these results have to be multiplied by n. 

V2(3+V3) J 
16(6-%/3) a 

for the square-kagome and 5 = S* = for the kagome 
lattice. 

We have considered only a special class of lattice de- 
formations (under which the independent localized one- 
magnon states survive). Although we are not able to 
prove rigorously that these lattice deformations are the 
most favorable, we have presented above nonrigorous ar- 
guments that these deformations take advantage of the 
localized magnons in an optimal way. However, we have 
rigorously shown that there exist lattice deformations 
which yield a gain in the total energy for large values 
of leading to a spin-Peierls instability of the lattice 
for an appropriate (large) magnetic field hi. 

To discuss the scenario of spin-Peierls instability more 
specifically we consider a magnetic field above the sat- 
uration field hi. For the corresponding fully polar- 



Minimal total energy is obtained for S 



S* 



ized ferromagnetic state a lattice distortion is not favor- 
able. Decreasing h till hi the homogeneous ferromag- 
netic state transforms into the "distorted magnon crys- 
tal" ; this transformation is accompanied by the afore- 
mentioned magnetization jump. On the basis of gen- 
eral arguments2i24 we expect that the "magnon crys- 
tal" state has gapped excitations and the system exhibits 
a magnetization plateau between hi and /12 < hi at 
= — rimax- To support this statement we cal- 
culate the plateau width Aft. — hi ^ h2 for finite systems 
of = 27,36,45,54 (kagome) and N = 24,30,48,54 
(square-kagome) for the undistorted lattice, where /12 is 



obtained by /12 = i?min(5'^ — 



) — E^in{S^ 



■^N ~ rimax ~ !)■ Using a finite-size extrapolation we 
find indeed evidence for a finite Ah = 0.07J {Ah = 
0.33J) for the kagome (square-kagome) lattice in the 
thermodynamic limit. This plateau width might be en- 
larged by distortions (see below). 

Now the question arises whether the lattice distortion 
under consideration is stable below this plateau, i.e., for 
< — rimax- We are not able to give a rigor- 
ous answer but can discuss the question again for fi- 
nite systems of size iV = 24, 30, 48, 54 (square-kagome) 
and N = 18,27,36,45,54 (kagome) with rimax distorted 
squares/hexagons. We calculate the magnetic energy for 
zero and small distortion parameter 5 for different values 
of S^- . Adopting for the magnetic energy the ansatz 



^^min(^^5) =^;min(5^0)+A(5P 



(9) 



and taking 5 of the order of 10"** we can estimate the 
exponent p from the numerical results. Evidently, the 
lattice may become unstable if p < 2 (and, of course, 
A <0) whereas p> 2 indicates lattice stability. 

Of course, the numerical results reproduce the ana- 
lytical findings reported above for 5*^ = ^A^ — fimax- 
More interesting is the sector of just below, i.e. 

= ^A^ — ?^max — 1- Remarkably both lattices be- 
have differently as h becomes smaller than /i2- For 
the finite square-kagome lattices considered we find p — 
1.001, 1.000, 1.000, 1.000 for A^ = 24,30,48,54, respec- 



tively, US' ^^N~ 



1. Moreover, p remains equal 



to 1 for smaller S'. Therefore, we conclude that the 
distorted square-kagome lattice remains stable for S"^ < 
— rimax- On the other hand, for finite kagome lat- 
tices we obtain p = 2.000, 1.002, 2.000, 1.998, 2.055 for 
A^= 18, 27, 36, 45, 54, respectively, ff 5^ = ^N~n,^^^-1. 
The "outrider" for N — 27 may be attributed to finite- 
size effects, indeed we have p = 1.994 for the next lower 
S' = — nmax — 2. Moreover, p remains about 2 for 
smaller 5^. We interpret small deviations from 2 also 
as finite-size effects and conclude that the spin-Peierls 
instability in the kagome lattice (within the adopted 
ansatz for the lattice deformation) is favorable only for 
iA^ — Tiinax < "S*^ < and the distortion disappears for 
h < /i2. 

The origin for the different behavior of the square- 
kagome and the kagome lattice below /12 can be at- 
tributed to the circumstance, that the square-kagome lat- 
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tice has non-equivalent NN bonds (namely bonds belong- 
ing to triangles and bonds belonging to squares) which 
carry different NN spin-spin correlations, but the kagome 
lattice has not. This difference in the NN spin-spin corre- 
lations leads to a special affinity of the magnetic system 
to the considered lattice distortions. 

Let us now briefly discuss the influence of the lattice 
distortion on the saturation field hi. Since the fully po- 
larized ferromagnetic state is not distorted, hi is shifted 

to higher values according to ''^^f =3-1- ^(^+'^) J. fQj. 

^ o J 32(6-\/3) a 

the square- kagome and ''^^ = 3 + ^ for the kagome 
lattice, i.e. the "distorted magnon crystal" remains sta- 
ble until hi{6*) > hi. On the other hand, starting at 
large magnetic fields h > hi the fully polarized ferromag- 
netic state remains (meta)stable until hi — 3J. Conse- 
quently, we have a hysteresis phenomenon in the vicinity 
of saturation field. 

We mention that our considerations basically remain 
unchanged for the anisotropic Hamiltonian with A ^ 
1. The presence of anisotropy leads only to quantita- 
tive changes in our results. For the sectors of be- 
low the localized magnon states we have checked that 
numerically for N — 18 (kagome) and N = 24 (square- 
kagome). For the sectors of with localized magnon 
states it becomes obvious from the change in magnetic 
energy of a localized magnon state due to distortions 
given by — (\/2 + (-\/3 — f ) A) SJ (square-kagome) 
and — (l + iA) SJ (kagome), where these expressions for 
A = f transform to the first terms in Eqs. (TJ and ||SJ|, 
respectively. 

From the experimental point of view the discussed ef- 
fect should most spectacularly manifest itself as a hys- 
teresis in the magnetization and the deformation of 
kagome-lattice antiferromagnets or the kagome- like mag- 
netic molecules in the vicinity of the saturation field. 

We wish to stress that the predicted spin-Peierls in- 
stability in high magnetic fields may appear in a whole 
class of frustrated quantum magnets in one, two and 
three dimensions hosting independent localized magnons 
provided it is possible to construct a lattice distortion 
preserving the symmetry of the localized-magnon cell. 



Moreover, the effect is not restricted to s = i and to 
isotropic Heisenberg systemsiiiSS. We mention, that we 
have checked explicitly that our results are valid for e.g. 
the kagome-like chain of Ref. or the dimer-plaquette 
chain (see Ref. and references therein). This fact 
certainly increases the chance to observe the predicted 
spin-Peierls instability. 

There is an increasing number of synthesized quan- 
tum frustrated kagome magnets^LSSiSSiSi. Though these 
available materials do not fit perfectly to an ideal kagome 
Heisenberg antiferromagnet the physical effects based on 
localized magnon states may survive in non-ideal geome- 
tries, see Ref. |^ Furthermore one needs compara- 
bly small exchange constants J to reach experimentally 
the saturation field where the structural instability oc- 
curs. A simple calculation leads to the relation /isat /Tesla 
^ 2.23 J/K for a spin-^ kagome Heisenberg antiferro- 
magnet. As a candidate for an experimental study at 
saturation field may serve the novel spin-| kagome like 
material Ba2Sn2ZnCr7pGaio-7p02^^ with a comparably 
small exchange constant of about J ~ 37 . . . 40K leading 
to a saturation field accessible with today available high- 
end experimental equipment. We may expect that fur- 
ther such materials will be synthesized and with our work 
we are pointing out that the efforts in this direction are 
worthwhile also because of a new effect: the spin-Peierls 
instability in a strong magnetic field. 

To summarize, we have reported a spin-Peierls in- 
stability in strong magnetic fields for several frustrated 
Heisenberg antiferromagnets hosting independent local- 
ized magnons. This spin-Peierls instability may or may 
not survive for smaller fields in dependence on details of 
lattice structure. In particular, for the Heisenberg an- 
tiferromagnet on the kagome lattice we have found ev- 
idence that the spin-Peierls instability breaking sponta- 
neously the translational symmetry of the kagome lattice 
appears only in a certain region of the magnetic field. 
The field dependence of the magnetization and the defor- 
mation in the vicinity of saturation displays a hysteresis. 

We thank A. Honecker, J. Schnack, S. L. Drechsler and 
D. Ihle for discussions and comments. The present study 
was supported by the DFG (project 436 UKR 17/17/03). 
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